Evolution of collision numbers for a chaotic gas dynamics 
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We put forward a conjecture of recurrence for a gas of hard spheres that collide elastically in a 
finite volume. The dynamics consists of a sequence of instantaneous binary collisions. We study how 
the numbers of coUisions of different pairs of particles grow as functions of time. We observe that 
these numbers can be represented as a time-integral of a function on the phase space. Assuming 
the results of the ergodic theory apply, we describe the evolution of the numbers by an effective 
Langevin dynamics. We use the facts that hold for these dynamics with probability one, in order 
to establish properties of a single trajectory of the system. We find that for any triplet of particles 
there will be an infinite sequence of moments of time, when the numbers of collisions of all three 
different pairs of the triplet will be equal. Moreover, any value of difference of collision numbers of 
pairs in the triplet will repeat indefinitely. On the other hand, for larger number of pairs there is 
but a finite number of repetitions. Thus the ergodic theory produces a limitation on the dynamics. 
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The systems of colliding hard-core particles take a spe- 
cial place in the theory of many-body systems. They were 
used extensively already by Boltzmann to study the fun- 
damental principles of the statistical physics These 
systems are unique as they allow insight into the basic 
properties of complex systems that are usually postu- 
lated. Namely, rigorous facts on the ergodicity 0, Q of 
these systems are known. The famous Boltzmann-Sinai 
hypothesis states that systems of an arbitrary number 

> 2 of elastic hard balls in a d— dimensional box with 
periodic boundary conditions (torus), d>2, are ergodic 
in the phase space region where the trivial conserved 
quantities of the system are constant [3]. The excep- 
tional feature of this hypothesis is that by today it can 
be considered as " almost proved" , in contrast to the er- 
godicity hypothesis on other chaotic systems where the 
proof is generally absent. The first rigorous result was 
obtained in 1970 by Sinai who proved that a system of 
two disks in a 2— dimensional torus is ergodic @. No- 
tably, this result shows amply that the thermodynamic 
limit of a large number of particles is not necessary for 
the ergodicity. The extension of the proof to an arbitrary 
number of particles and arbitrary d is " almost" complete 
by now, see e. g. @ and references therein, and also 
0, Q . As the system is equivalent to a billiard 3 " ^ 
single particle colliding elastically against the boundary 
of a certain manifold - we will call it below a " billiard" . 

In this Rapid Communication we use the ergodic the- 
ory not to deal with the equilibrium properties of the 
billiard, but rather to extract information on the struc- 
ture of dynamics at finite times (this means considering 
non-equilibrium phenomena, as the equilibrium statisti- 
cal physics, within the approach of the ergodic theory, 
describes the infinite-time averages). The dynamics are 
a sequence of the events of binary collisions of particles. 
This sequence is an ordered list of pairs of particles, say 
(1,2) (5,6) (7,8)... meaning that first particles 1 and 2 
collided, then 5 and 6, then 7 and 8 and so on. Due to 
chaos this list looks like a random sequence of pairs. Here 



we find a deterministic constraint on that sequence. 

The sequence of collisions Ci is a discrete process tak- 
ing K{K — l)/2 values corresponding to the pairs of the 
system with K particles. This process is not a Bernoulli 
scheme, i. e. Ci and Q+i, determined by Newton's law, 
are not independent. However, the Bernoulli property, 
which is likely to hold for the considered system [l^, llll , 
would suggest that Ci and Ci+m do become independent 
in the limit m — )■ oo. Thus a coarse-graining of c; over a 
sufficient number of steps would produce a process which 
consecutive steps are already independent to a good ap- 
proximation, producing a random walk. Below we per- 
form the analysis of the coarse-grained sequence in the 
formulation that we found convenient. 

We consider the numbers of collisions of pairs of par- 
ticles up to a time t, cf [12 • Roughly, for sufficiently 
large t one can split the considered time-interval into 
sub-intervals with approximately independent numbers 
of collisions, so the total number of collisions is a sum 
of many i. i. d. random variables and one can use the 
central limit theorem (CLT). Below we formulate this as- 
sumption using the ergodic theory. Then, the numbers 
of collisions of different pairs can be described using an 
effective system of Langevin equations. The sequence of 
collision events following from the Newtonian dynamics 
is statistically indistinguishable from a realization of the 
random Langevin dynamics, cf. [lo| . Here the statistics 
is defined by the volume of the initial conditions in the 
phase space corresponding to the considered property. 

The effective Langevin description allows to put for- 
ward a conjecture of recurrence that seems to be new. 
For any triplet of particles in the billiard, there is an in- 
finite sequence of times at which the collision numbers 
of the corresponding three pairs of particles are equal. 
For AT = 3 this conjecture covers all the particles of the 
system. Furthermore, any prescribed value of two dif- 
ferences of collisions will repeat infinitely many times. 
In contrast, a prescribed value of three and more differ- 
ences, existing for A' > 3, repeats but a finite number of 
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times. Significantly, the statements are dynamical and 
they hold for almost every trajectory of the system (i. 
e. with a possible exception of a set of trajectories the 
volume of the initial conditions of which is zero). 

The analysis below applies to billiards both in I? = 2 
and D = 3 cases, where D is the space dimension, with 
ramifications in D = 2 case entailed by the slow, non- 
integrable decay of the correlation tails [l^, [l^ . For def- 
initeness one can think of hard balls of diameter d that 
collide in a square cube with periodic boundary condi- 
tions. The dynamics is a succession of the discrete events 
of binary collisions. Starting from a given initial condi- 
tion for particles' positions and velocities, one determines 
the time to the next collision and the pair that is going 
to collide. Pushing then the particles' positions and the 
velocities to the time after the collision, one iterates the 
procedure. We designate the number of collisions of par- 
ticles i and j that occurred up to time t by Nij (t) . The 
important observation at the basis of the analysis be- 
low is that Nij (i) can be represented as an integral of a 
function on the phase space. To provide this representa- 
tion we first consider a formal representation of Nij (t) in 
terms of the Dirac i5— function. We designate the collision 



The form ([2]) of Nij (t) suggests the use of the results of 
the ergodic theory. Ergodicity guarantees the existence 
of the average collision rate limt _^oo Nfj [t)/t= v'^ , 
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where the angular brackets stand for the average over the 
microcanonical ensemble, and z^' does not depend on ij . 
It seems safe to assume that v'^ has a finite limit v as 
e — ?► 0. While the above applies to a single trajectory, 
here we address the behavior of the trajectories statis- 
tically. The statistics is defined by picking the initial 
condition at random in the allowed region of the phase 
space. The latter is defined by the trivial conserved quan- 
tities of the system (for torus energy and momentum), i. 
e. we consider the microcanonical ensemble. 

One expects that correlations of at large times de- 
cay as i"^/^, where D is the space dimension [l^ [l^ . 
For D ~ 3 the integral of t^^/^ converges and one 
may assume ^fj (t) effectively has a finite correlation time 
T^^j. < cx), cf. e. g. 0- Here we make the main assump- 
tion underlying the conjecture proposed in this work. We 



times of the pair i and j by t^^ , and note that {t) - (f is assume the correlation time t^^^ has a finite limit r< 



a non- negative function of t that vanishes only at t = tfj . 
Introducing F{r, v) = 2S [r^ — cP] \r ■ v\, write formally 
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where {t) and Vij (t) are the relative distance and ve- 
locity of particles i and j, respectively. Using the above, 
one may introduce a rigorous representation of Nijit) 
by any family of smooth functions used to represent the 
(5— function. For example, using exponential profiles, and 
the fact that r| {t) — d^ > except at t 
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e — ?► 0. While this assumption seems very plausible (the 
representation seems to work for the dilute gas [l3l where 
no increase of the correlation time is known to us), it is 
this assumption that allows us to circumvent the singu- 
larity in the representation ([IJ. Then at i 3> Tcor the 
numbers N^j{t) are sums of roughly t/rco,- ^ 1 indepen- 
dent random variables and one can use the CLT, giving 
the following Gaussian approximation to the probability 
density function (PDF) P{{Nij},t) of Nij{t) at t > Tcor- 
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exp [-{N, - t.Or,r.;^^(jV,„„ _ „t)/4t] 
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where [t) are the regularized collision rates, 

it^{t) ^ 2e-V,, • exp [-e'^ {r^ - d')] , (3) 

rigorously. To see how the above works, single out a 
neighborhood of the collision {tfj — S, tfj +6), where S is so 
small one can assume the particles i and j do not collide 
with other particles during that interval. The quantity S 
depends on properties of the system, like concentration. 
Considering a sufficiently small e <^ S we have 



where the summation over repeated indices is assumed. 
The dispersion matrix Tij^mn describes the fluctuations 
of the collision rates £,ij (t) and it is given by (the double 
brackets stand for the dispersion so that for any random 
variables x, y we have {{xy)) = (xy) — {x){y)): 
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where we used that • Vij < before the collision and 
Tij ■ Vij > after the collision. Note the latter representa- 
tion can be generalized straightforwardly to a dilute gas 
of particles with short-range interactions [17j . 



The use of the CLT above neglects the tails in 
P{{Nij},t). For certain billiards these tails can be even 
algebraic, see e. g. {l3|. This is not a limitation for the 
analysis below, that concerns the probability of fixed val- 
ues of Nij{t) — N„m{t) at large t. At a sufficiently large 
t, any given value belongs to the bulk of the PDF of the 
differences and is describable by the CLT. 

The simplest system for which the above relations ap- 
ply is the D = 3 system of 2 balls in a torus. In this case 
we have but one pair of particles, so that the number of 
collisions N that occurred by the time t obeys 
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The above relation appears to be a fundamental result 
on a very basic system and it demands further studies, 
both theoretical and numerical that we postpone for fur- 
ther work We note the same statistical distribution 
P{{Nij},t) would result for the Langevin dynamics 

^ = QAt), (C.,) = ^, (7) 
((C^J WC™«(t'))) = 2ry- ,„„,5(i - 1'). (8) 

Thus the above stochastic dynamics gives the effective 
description of the collision numbers Nij{t) in quite the 
same sense as the usual Langevin dynamics does: if one 
considers the dynamics over the temporal scale of coarse- 
graining that is much larger than Tcor then the two dy- 
namics are statistically equivalent. 

Eqs. (0)-® differ from similar results of the equilib- 
rium statistical physics. There the Langevin dynamics 
describes macroscopic quantities determined by a large 
number of particles. Here the quantities Nijit) are not 
macroscopic since the result holds even for the system of 
two particles (at least for the torus) . The macroscopic na- 
ture of the law is due to the consideration of the dynam- 
ics on a large time-scale. While Nij{t) depend strongly 
on the details of the Newtonian mechanics at time-scales 
< Tcori on a larger time-scale the dynamics forgets the 
details of the mechanism of collisions and Nij{t) are ef- 
fectively Brownian motions with non-zero mean. 

Using the above Langevin dynamics and facts that hold 
for the Brownian motion with probability one, we make 
deterministic predictions on the billiard. We note that 
the average rates of growth of TV^ {t) are equal, so the dif- 
ferences Nij = Nij — Ni2 are regular Brownian motions, 

— ^ — LOij, {iOij{t)uJmn{t')) = 2Dij^mnS{t — t'), 
Dij^mn = ^ij.mn + 1^12,12 ^ ^ rmri42- (9) 

The non-diagonality of Dij^mn is not important for our 
considerations below. We note however that using that 
Dij^mn is symmetric it is always possible to pass to ro- 
tated Nij that perform independent Brownian motions. 

We use the familiar fact that Brownian motion returns 
to the origin with unit probability for dimension lower or 
equal to the critical dimension 2. In D > 2 the return is 
probabilistic - there is a finite probability of return which 
is strictly less than one [l^, so the number of returns to 
the origin is always finite. In the light of this, the differ- 
ences Nij, that exist for systems with K > 3, are seen to 
be special in the case K = 3. Here the Brownian motion 
(iVi3, A^23) is two-dimensional. We arrive at the following 
conjecture of recurrence: for billiards with three particles 
there will almost always (i. e. with a possible exception 
of trajectories which initial conditions have zero volume) 
be a time when the numbers of collisions of all three pairs 
will equalize, Ni2{t) = N23{t) = Ni2{t). Moreover, the 
system will be getting back to these equalized states an 
infinite number of times. Here it should be clear that 



the equality sign should be understood with a finite ac- 
curacy following both from the approximate nature of 
the Langevin equation and from the fact that the Brow- 
nian motion in = 2 is only neighborhood-recurrent 
and not point- recurrent as in = 1 (l6j . This limi- 
tation is not important qualitatively, since Nij{t) grow 
with time indefinitely, so the finite accuracy is irrele- 
vant at large times. The result is quite distinct from 
the familiar Poincare recurrence theorem where the sys- 
tem gets back to the neighborhood of the same point in 
the phase space: the two recurrences are generally unre- 
lated. Furthermore, based on the fact that the Brownian 
motion visits neighborhood of every point in the plane 
an indefinite number of times, the conjecture can be ex- 
tended to the statement that every possible combination 
of [Nisit) - Ni2{t), N23{t) - Ni2{t)] is going to occur and 
then recur an infinite number of times. 

The above conjecture of recurrence does not hold for 
systems with K > 4. Here with a finite probability after 
t = there will never be again a situation where all Nij (t) 
are equal. The volume fraction of the initial conditions 
in the phase space for which all Nij (t) get equal at some 
time t ^ Tcor is strictly less than one for if > 4 and it de- 
creases as K grows (for the explicit formula see [11] )• A 
limited version of the recurrence hypothesis holds, stat- 
ing that for any three variables Nij , N^i and Nmn there 
is an infinite sequence of times at which any two dif- 
ferences of these numbers will take a fixed, preassigned 
value. However, with probability one, three and more 
linearly independent differences of collision numbers, de- 
finable for a number of particles larger than three, will 
repeat but a finite number of times. Finally, the return 
property in I? = 1, gives that the sequence of times at 
which Nijit) = Nmn(t) is infinite for any z, j, m and n. 

We now rederive the conjecture qualitatively. We con- 
sider particles 1, 2 and 3 and study the D = 2 vector 
(A^i3, N23)- At large t the probability that iVia equals to 
zero decays as t~^/^. Assuming effective independence of 
Nisit) and N23{t) we conclude that the probability that 
both latter functions equal zero decays as t~^. Since the 
latter is non-integrable, the probability that there exists 
a finite time t such that Ni3{t) — N23{t) = is one. 
For three functions Nij{t) the corresponding probability 
would decay as t~'^/^ and the recurrence probability is 
strictly less than one. This argument can be generalized 
to the recurrence of (iVi3,iV23) = (m,n) for any m, n. 

The derivation can be generalized to D = 2, where the 
usual Langevin dynamics does not hold. Such a general- 
ization is important because billiards on D = 2 tori were 
studied extensively. For D = 2 the correlations of ^ij are 
expected [H, to decay as I/t making Tij^mn in Eq- © 
divergent. As a result the CLT cannot be used, how- 
ever the recurrence still holds. The dispersion of Nij{t) 
grows at tint, so the probability that Nij{t) = de- 
cays as (tlni)~^/^. The logarithmic correction does not 
change the convergence/divergence of the corresponding 
integrals and the argument above can be repeated. Thus 
we expect the conjecture also to hold in D = 2. 
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We showed a use of the ergodic theory to approach a 
single reahzation of a chaotic dynamics. The basic obser- 
vation is that under the assumption of a finite correlation 
time, natural for the considered systems, the PDF of the 
numbers of collisions of pairs of particles is described by 
the CLT. This allows to introduce Langevin's equations 
that provide an effective description of the dynamics of 
those numbers. Based on this description, the behav- 
ior of the realizations of which is well-known, one can 
draw conclusions on the behavior of the billiard's trajec- 
tories. For this one concentrates on facts that hold for 
the Langevin dynamics with probability one. Using the 
return properties of the random walks, we demonstrated 
that for every triplet of particles in the billiard there 
is an infinite sequence of times where the differences of 
the collision numbers of all three pairs take an arbitrary 
fixed value. Furthermore, any value of three and more 
linearly independent differences of collision numbers, de- 
finable for K > 3, repeats but a finite number of times. 

The conjecture involves the continuous approximation 
to the discrete process Ci via the approximate dynam- 
ics of Nij{t). It would be advantageous to study the 
discrete random walk Ci directly. This would demand 
the consideration of the one-dimensional random walk 
-^jj(^) ~ Nmn{t) and the two-dimensional discrete ran- 
dom walk [Nij (t) - N,nn{t),Npr{t) - Nki{t)]. Both walks 
occur on the lattice with step size one. However, there 
are two major differences from the usual discrete random 
walk. First, the times between different steps of the walk 
are random. Second, both, the time intervals between 
one step of the walk and the next one, and Ci and c^+i are 
not independent. Thus the description of the appropriate 
random walk would demand either the consideration of 
detailed joint PDFs of time-intervals between the steps, 
or introducing a certain coarse-graining in time making 
the subsequent steps of the walk independent. The lat- 
ter approach would involve a continuous approximation 
to Nij{t) and it can be said that our work does just that. 



The described phenomena apply to any number of par- 
ticles larger than one as long as the ergodic theory ap- 
plies. Since they concern large numbers of collisions of 
the same pairs of particles, their practical observation 
(either numerical or experimental) demands considering 
systems with a relatively small number of particles. The 
results might also apply to the Lyapunov modes of large 
systems, that involve a small number of particles [l^ . 

The hypothesis applies to almost every trajectory only. 
In particular, we do not consider those initial conditions 
for which the hard-sphere dynamics is ill-defined, which 
phase space volume was shown to vanish (Toj . 

Within the approximation used in the Boltzmann 
equation, the collisions in the gas of particles with short- 
range interactions are considered instantaneous. Thus 
the conjecture proposed here can be expected to hold for 
any gas of particles, where similar statistical properties of 
Nij{t) hold. In particular, our derivation of the effective 
Langevin description, allows direct generalization to the 
gas [for the gas, one can again introduce the time-scale 5 
where the colliding pair is "isolated" and the appropriate 
spatial scale We will report elsewhere the numeri- 
cal studies of the conjecture proposed here, both for hard 
spheres and for the gas with short-range interactions [13] . 

In this work we used the ergodic theory to derive quite 
detailed dynamical properties of almost all trajectories of 
the billiard. We showed that any prescribed value of two 
differences of collision numbers will repeat indefinitely. 
On the other hand, any value of three and more differ- 
ences of collision numbers repeats but a finite number 
of times. These properties seem to be non-obvious from 
dynamical standpoint. The study of the underlying as- 
sumption of the asymptotic independence of Ci and Ci+m 
at large to, both numerical and theoretical, is the subject 
for further analysis [l7|. 

We thank N. Chernov and N. Simanyi for very useful 
remarks that helped to improve the paper significantly. 
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